An Endosomal Acid-Regulatory Feedback System Rewires Cytosolic cAMP Metabolism and Drives Tumor Progression.

Although suppressed cAMP levels have been linked to cancer for nearly five decades, the molecular basis remains uncertain. Here, we identify endosomal pH as a novel regulator of cytosolic cAMP homeostasis and a promoter of transformed phenotypic traits in colorectal cancer. Combining experiments and computational analysis, we show that the Na+/H+ exchanger NHE9 contributes to proton leak and causes luminal alkalinization, which induces resting [Ca2+], and in consequence, represses cAMP levels, creating a feedback loop that echoes nutrient deprivation or hypoxia. Higher NHE9 expression in cancer epithelia is associated with a hybrid epithelial-mesenchymal (E/M) state, poor prognosis, tumor budding, and invasive growth in vitro and in vivo. These findings point to NHE9-mediated cAMP suppression as a pseudostarvation-induced invasion state and potential therapeutic vulnerability in colorectal cancer. Our observations lay the groundwork for future research into the complexities of endosome-driven metabolic reprogramming and phenotype switching and the biology of cancer progression. IMPLICATIONS: Endosomal pH regulator NHE9 actively controls cytosolic Ca2+ levels to downregulate the adenylate cyclase-cAMP system, enabling colorectal cancer cells to acquire hybrid E/M characteristics and promoting metastatic progression.

Multi-modal transcriptomic analysis unravels enrichment of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer.

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.

Aberrations in ion channels interacting with lipid metabolism and epithelial-mesenchymal transition in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is the most prevalent malignant gastrointestinal tumor. Ion channels contribute to tumor growth and progression through interactions with their neighboring molecules including lipids. The dysregulation of membrane ion channels and lipid metabolism may contribute to the epithelial-mesenchymal transition (EMT), leading to metastatic progression. Herein, transcriptome profiles of patients with ESCC were analyzed by performing differential gene expression and weighted gene co-expression network analysis to identify the altered ion channels, lipid metabolism- and EMT-related genes in ESCC. A total of 1,081 differentially expressed genes, including 113 ion channels, 487 lipid metabolism-related, and 537 EMT-related genes, were identified in patients with ESCC. Thereafter, EMT scores were correlated with altered co-expressed genes. The altered co-expressed genes indicated a correlation with EMT signatures. Interactions among 22 ion channels with 3 hub lipid metabolism- and 13 hub EMT-related proteins were determined using protein-protein interaction networks. A pathway map was generated to depict deregulated signaling pathways including insulin resistance and the estrogen receptor-Ca2+ signaling pathway in ESCC. The relationship between potential ion channels and 5-year survival rates in ESCC was determined using Kaplan-Meier plots and Cox proportional hazard regression analysis. Inositol 1,4,5-trisphosphate receptor type 3 (ITPR3) was found to be associated with poor prognosis of patients with ESCC. Additionally, drugs interacting with potential ion channels, including GJA1 and ITPR3, were identified. Understanding alterations in ion channels with lipid metabolism and EMT in ESCC pathophysiology would most likely provide potential targets for the better treatment of patients with ESCC.

Unique sphingolipid signature identifies luminal and triple-negative breast cancer subtypes.

Breast cancer (luminal and triple-negative breast cancer [TNBC]) is the most common cancer among women in India and worldwide. Altered sphingolipid levels have emerged as a common phenomenon during cancer progression. However, these alterations are yet to be translated into robust diagnostic and prognostic markers for cancer. Here, we present the quantified sphingolipids of tumor and adjacent-normal tissues from patients of luminal (n = 70) and TNBC (n = 42) subtype from an Indian cohort using targeted liquid chromatography mass spectrometry. We recorded unique sphingolipid profiles that distinguished luminal and TNBC tumors in comparison to adjacent normal tissue by six-sphingolipid signatures. Moreover, systematic comparison of the profiles of luminal and TNBC tumors provided a unique five-sphingolipid signature distinguishing the two subtypes. We further identified key sphingolipids that can stratify grade II and grade III tumors of luminal and TNBC subtype as well as their lymphovascular invasion status. Therefore, we provide the right evidence to develop these candidate sphingolipids as widely acceptable marker/s capable of diagnosing luminal vs TNBC subtype of breast cancer, and predicting the disease severity by identifying the tumor grade.

PD-L1 Activity Is Associated with Partial EMT and Metabolic Reprogramming in Carcinomas.

Immune evasion and metabolic reprogramming are hallmarks of cancer progression often associated with a poor prognosis and frequently present significant challenges for cancer therapies. Recent studies have highlighted the dynamic interaction between immunosuppression and the dysregulation of energy metabolism in modulating the tumor microenvironment to promote cancer aggressiveness. However, a pan-cancer association among these two hallmarks, and a potent common driver for them-epithelial-mesenchymal transition (EMT)-remains to be done. This meta-analysis across 184 publicly available transcriptomic datasets as well as The Cancer Genome Atlas (TCGA) data reveals that an enhanced PD-L1 activity signature along with other immune checkpoint markers correlate positively with a partial EMT and an elevated glycolysis signature but a reduced OXPHOS signature in many carcinomas. These trends were also recapitulated in single-cell, RNA-seq, time-course EMT induction data across cell lines. Furthermore, across multiple cancer types, concurrent enrichment of glycolysis and PD-L1 results in worse outcomes in terms of overall survival as compared to enrichment for only PD-L1 activity or expression. These results highlight potential functional synergy among these interconnected axes of cellular plasticity in enabling metastasis and multi-drug resistance in cancer.

Analysis of Yes-Associated Protein-1 (YAP1) Target Gene Signature to Predict Progressive Breast Cancer.

Breast cancers are treated according to the ER/PR or HER2 expression and show better survival outcomes with targeted therapy. Triple-negative breast cancers (TNBCs) with a lack of expression of ER/PR and HER2 are treated with systemic therapy with unpredictable responses and outcomes. It is essential to investigate novel markers to identify targeted therapies for TNBC. One such marker is YAP1, a transcription co-activator protein that shows association with poor prognosis of breast cancer. YAP1 transcriptionally regulates the expression of genes that drive the oncogenic phenotypes. Here, we assess a potential YAP target gene signature to predict a progressive subset of breast tumors from METABRIC and TCGA datasets. YAP1 target genes were shortlisted based on expression correlation and concordance with YAP1 expression and significant association with survival outcomes of patients. Hierarchical clustering was performed for the shortlisted genes. The utility of the clustered genes was assessed by survival analysis to identify a recurring subset. Expression of the shortlisted target genes showed significant association with survival outcomes of HER2-positive and TNBC subset in both datasets. The shortlisted genes were verified using an independent dataset. Further validation using IHC can prove the utility of this potential prognostic signature to identify a recurrent subset of HER2-positive and TNBC subtypes.

In Silico Analysis of Ion Channels and Their Correlation with Epithelial to Mesenchymal Transition in Breast Cancer.

Uncontrolled growth of breast cells due to altered gene expression is a key feature of breast cancer. Alterations in the expression of ion channels lead to variations in cellular activities, thus contributing to attributes of cancer hallmarks. Changes in the expression levels of ion channels were observed as a consequence of EMT. Additionally, ion channels were reported in the activation of EMT and maintenance of a mesenchymal phenotype. Here, to identify altered ion channels in breast cancer patients, differential gene expression and weighted gene co-expression network analyses were performed using transcriptomic data. Protein-protein interactions network analysis was carried out to determine the ion channels interacting with hub EMT-related genes in breast cancer. Thirty-two ion channels were found interacting with twenty-six hub EMT-related genes. The identified ion channels were further correlated with EMT scores, indicating mesenchymal phenotype. Further, the pathway map was generated to represent a snapshot of deregulated cellular processes by altered ion channels and EMT-related genes. Kaplan-Meier five-year survival analysis and Cox regressions indicated the expression of CACNA1B, ANO6, TRPV3, VDAC1 and VDAC2 to be potentially associated with poor survival. Deregulated ion channels correlate with EMT-related genes and have a crucial role in breast cancer-associated tumorigenesis. Most likely, they are potential candidates for the determination of prognosis in patients with breast cancer.

Quantifying the Patterns of Metabolic Plasticity and Heterogeneity along the Epithelial-Hybrid-Mesenchymal Spectrum in Cancer.

Cancer metastasis is the leading cause of cancer-related mortality and the process of the epithelial-to-mesenchymal transition (EMT) is crucial for cancer metastasis. Both partial and complete EMT have been reported to influence the metabolic plasticity of cancer cells in terms of switching among the oxidative phosphorylation, fatty acid oxidation and glycolysis pathways. However, a comprehensive analysis of these major metabolic pathways and their associations with EMT across different cancers is lacking. Here, we analyse more than 180 cancer cell datasets and show the diverse associations of these metabolic pathways with the EMT status of cancer cells. Our bulk data analysis shows that EMT generally positively correlates with glycolysis but negatively with oxidative phosphorylation and fatty acid metabolism. These correlations are also consistent at the level of their molecular master regulators, namely AMPK and HIF1α. Yet, these associations are shown to not be universal. The analysis of single-cell data for EMT induction shows dynamic changes along the different axes of metabolic pathways, consistent with general trends seen in bulk samples. Further, assessing the association of EMT and metabolic activity with patient survival shows that a higher extent of EMT and glycolysis predicts a worse prognosis in many cancers. Together, our results reveal the underlying patterns of metabolic plasticity and heterogeneity as cancer cells traverse through the epithelial-hybrid-mesenchymal spectrum of states.

Semicoordinated allelic-bursting shape dynamic random monoallelic expression in pregastrulation embryos.

Recently, allele-specific single-cell RNA-seq analysis has demonstrated widespread dynamic random monoallelic expression of autosomal genes (aRME) in different cell types. However, the prevalence of dynamic aRME during pregastrulation remains unknown. Here, we show that dynamic aRME is widespread in different lineages of pregastrulation embryos. Additionally, the origin of dynamic aRME remains elusive. It is believed that independent transcriptional bursting from each allele leads to dynamic aRME. Here, we show that allelic burst is not perfectly independent; instead it happens in a semicoordinated fashion. Importantly, we show that semicoordinated allelic bursting of genes, particularly with low burst frequency, leads to frequent asynchronous allelic bursting, thereby contributing to dynamic aRME. Furthermore, we found that coordination of allelic bursting is lineage specific and genes regulating the development have a higher degree of coordination. Altogether, our study provides significant insights into the prevalence and origin of dynamic aRME and their developmental relevance during early development.

An Integrative Systems Biology Approach Identifies Molecular Signatures Associated with Gallbladder Cancer Pathogenesis.

Gallbladder cancer (GBC) has a lower incidence rate among the population relative to other cancer types but is a major contributor to the total number of biliary tract system cancer cases. GBC is distinguished from other malignancies by its high mortality, marked geographical variation and poor prognosis. To date no systemic targeted therapy is available for GBC. The main objective of this study is to determine the molecular signatures correlated with GBC development using integrative systems level approaches. We performed analysis of publicly available transcriptomic data to identify differentially regulated genes and pathways. Differential co-expression network analysis and transcriptional regulatory network analysis was performed to identify hub genes and hub transcription factors (TFs) associated with GBC pathogenesis and progression. Subsequently, we assessed the epithelial-mesenchymal transition (EMT) status of the hub genes using a combination of three scoring methods. The identified hub genes including, CDC6, MAPK15, CCNB2, BIRC7, L3MBTL1 were found to be regulators of cell cycle components which suggested their potential role in GBC pathogenesis and progression.

Epithelial-to-Mesenchymal Transition Enhances Cancer Cell Sensitivity to Cytotoxic Effects of Cold Atmospheric Plasmas in Breast and Bladder Cancer Systems.

Cold atmospheric plasma (CAP) has emerged as a highly selective anticancer agent, most recently in the form of plasma-activated medium (PAM). Since epithelial-mesenchymal transition (EMT) has been implicated in resistance to various cancer therapies, we assessed whether EMT status is associated with PAM response. Mesenchymal breast cancer cell lines, as well as the mesenchymal variant in an isogenic EMT/MET human breast cancer cell system (PMC42-ET/LA), were more sensitive to PAM treatment than their epithelial counterparts, contrary to their responses to other therapies. The same trend was seen in luminal muscle-invasive bladder cancer model (TSU-Pr1/B1/B2) and the non-muscle-invasive basal 5637 bladder cancer cell line. Three-dimensional spheroid cultures of the bladder cancer cell lines were less sensitive to the PAM treatment compared to their two-dimensional counterparts; however, incrementally better responses were again seen in more mesenchymally-shifted cell lines. This study provides evidence that PAM preferentially inhibits mesenchymally-shifted carcinoma cells, which have been associated with resistance to other therapies. Thus, PAM may represent a novel treatment that can selectively inhibit triple-negative breast cancers and a subset of aggressive bladder cancers, which tend to be more mesenchymal. Our approach may potentially be utilized for other aggressive cancers exhibiting EMT and opens new opportunities for CAP and PAM as a promising new onco-therapy.

A mechanistic model captures the emergence and implications of non-genetic heterogeneity and reversible drug resistance in ER+ breast cancer cells.

Resistance to anti-estrogen therapy is an unsolved clinical challenge in successfully treating ER+ breast cancer patients. Recent studies have demonstrated the role of non-genetic (i.e. phenotypic) adaptations in tolerating drug treatments; however, the mechanisms and dynamics of such non-genetic adaptation remain elusive. Here, we investigate coupled dynamics of epithelial-mesenchymal transition (EMT) in breast cancer cells and emergence of reversible drug resistance. Our mechanism-based model for underlying regulatory network reveals that these two axes can drive one another, thus enabling non-genetic heterogeneity in a cell population by allowing for six co-existing phenotypes: epithelial-sensitive, mesenchymal-resistant, hybrid E/M-sensitive, hybrid E/M-resistant, mesenchymal-sensitive and epithelial-resistant, with the first two ones being most dominant. Next, in a population dynamics framework, we exemplify the implications of phenotypic plasticity (both drug-induced and intrinsic stochastic switching) and/or non-genetic heterogeneity in promoting population survival in a mixture of sensitive and resistant cells, even in the absence of any cell-cell cooperation. Finally, we propose the potential therapeutic use of mesenchymal-epithelial transition inducers besides canonical anti-estrogen therapy to limit the emergence of reversible drug resistance. Our results offer mechanistic insights into empirical observations on EMT and drug resistance and illustrate how such dynamical insights can be exploited for better therapeutic designs.

Transcriptomic-Based Quantification of the Epithelial-Hybrid-Mesenchymal Spectrum across Biological Contexts.

Epithelial-mesenchymal plasticity (EMP) underlies embryonic development, wound healing, and cancer metastasis and fibrosis. Cancer cells exhibiting EMP often have more aggressive behavior, characterized by drug resistance, and tumor-initiating and immuno-evasive traits. Thus, the EMP status of cancer cells can be a critical indicator of patient prognosis. Here, we compare three distinct transcriptomic-based metrics-each derived using a different gene list and algorithm-that quantify the EMP spectrum. Our results for over 80 cancer-related RNA-seq datasets reveal a high degree of concordance among these metrics in quantifying the extent of EMP. Moreover, each metric, despite being trained on cancer expression profiles, recapitulates the expected changes in EMP scores for non-cancer contexts such as lung fibrosis and cellular reprogramming into induced pluripotent stem cells. Thus, we offer a scoring platform to quantify the extent of EMP in vitro and in vivo for diverse biological applications including cancer.

Immunosuppressive Traits of the Hybrid Epithelial/Mesenchymal Phenotype.

Recent preclinical and clinical data suggests enhanced metastatic fitness of hybrid epithelial/mesenchymal (E/M) phenotypes, but mechanistic details regarding their survival strategies during metastasis remain unclear. Here, we investigate immune-evasive strategies of hybrid E/M states. We construct and simulate the dynamics of a minimalistic regulatory network encompassing the known associations among regulators of EMT (epithelial-mesenchymal transition) and PD-L1, an established immune-suppressor. Our simulations for the network consisting of SLUG, ZEB1, miR-200, CDH1 and PD-L1, integrated with single-cell and bulk RNA-seq data analysis, elucidate that hybrid E/M cells can have high levels of PD-L1, similar to those seen in cells with a full EMT phenotype, thus obviating the need for cancer cells to undergo a full EMT to be immune-evasive. Specifically, in breast cancer, we show the co-existence of hybrid E/M phenotypes, enhanced resistance to anti-estrogen therapy and increased PD-L1 levels. Our results underscore how the emergent dynamics of interconnected regulatory networks can coordinate different axes of cellular fitness during metastasis.

Single-Cell Analysis Reveals Partial Reactivation of X Chromosome instead of Chromosome-wide Dampening in Naive Human Pluripotent Stem Cells.

Recently, a unique form of X chromosome dosage compensation has been demonstrated in human preimplantation embryos, which happens through the dampening of X-linked gene expression from both X chromosomes. Subsequently, X chromosome dampening has also been demonstrated in female human pluripotent stem cells (hPSCs) during the transition from primed to naive state. However, the existence of dampened X chromosomes in both embryos and hPSCs remains controversial. Specifically, in preimplantation embryos it has been shown that there is inactivation of X chromosome instead of dampening. Here, we performed allelic analysis of X-linked genes at the single-cell level in hPSCs and found that there is partial reactivation of the inactive X chromosome instead of chromosome-wide dampening upon conversion from primed to naive state. In addition, our analysis suggests that the reduced X-linked gene expression in naive hPSCs might be the consequence of erasure of active X chromosome upregulation.

Novel simulated moving-bed adsorber for the fractionation of gas mixtures.

The separation of propylene-propane mixture is an energy intensive operation commercially practiced using cryogenic distillation. The separation by pressure swing adsorption has been studied as an alternative. A fixed-bed pressure swing adsorption yields the heavy component as a pure product. The product recovery and the productivity are not high. In a moving-bed process, because of the counter-current solid-gas contact, the separation achieved is similar to that of the fractionation by distillation. Although the moving-bed operation offers the upper limit for the performance of a cyclic adsorptive process, due to mechanical complexities in the handling of solids the 'simulated' moving-bed is preferred. By moving the inlet and outlet ports of streams located along the length of the bed, a moving-bed process can be realized in a fixed bed. We describe here a 'moving-port' system which permits injection or withdrawal of the fluid along the axial direction in a fixed bed. A fixed bed embedded with the moving-port systems emulates a simulated moving-bed adsorber. The proposed adsorber can fractionate a binary gas mixture into two product streams with high purities. It is similar to the Sorbex process of UOP but does not have the eluent as an additional separating agent. A parametric study indicates that high purity products and a higher productivity by an order of magnitude can be achieved with simulated moving-beds compared to the fixed beds.